RrA, an enzyme from Rhodospirillum rubrum, is a prototype of a new family of short-chain L-asparaginases.

L-Asparaginases (ASNases) catalyze the hydrolysis of L-Asn to L-Asp and ammonia. Members of the ASNase family are used as drugs in the treatment of leukemia, as well as in the food industry. The protomers of bacterial ASNases typically contain 300-400 amino acids (typical class 1 ASNases). In contrast, the chain of ASNase from Rhodospirillum rubrum, reported here and referred to as RrA, consists of only 172 amino acid residues. RrA is homologous to the N-terminal domain of typical bacterial class 1 ASNases and exhibits millimolar affinity for L-Asn. In this study, we demonstrate that RrA belongs to a unique family of cytoplasmic, short-chain ASNases (scASNases). These proteins occupy a distinct region in the sequence space, separate from the regions typically assigned to class 1 ASNases. The scASNases are present in approximately 7% of eubacterial species, spanning diverse bacterial lineages. They seem to be significantly enriched in species that encode for more than one class 1 ASNase. Here, we report biochemical, biophysical, and structural properties of RrA, a member of scASNases family. Crystal structures of the wild-type RrA, both with and without bound L-Asp, as well as structures of several RrA mutants, reveal topologically unique tetramers. Moreover, the active site of one protomer is complemented by two residues (Tyr21 and Asn26) from another protomer. Upon closer inspection, these findings clearly outline scASNases as a stand-alone subfamily of ASNases that can catalyze the hydrolysis of L-Asn to L-Asp despite the lack of the C-terminal domain that is present in all ASNases described structurally to date.

Competition between electrostatic interactions and halogen bonding in the protein-ligand system: structural and thermodynamic studies of 5,6-dibromobenzotriazole-hCK2α complexes.

CK2 is a member of the CMGC group of eukaryotic protein kinases and a cancer drug target. It can be efficiently inhibited by halogenated benzotriazoles and benzimidazoles. Depending on the scaffold, substitution pattern, and pH, these compounds are either neutral or anionic. Their binding poses are dictated by a hydrophobic effect (desolvation) and a tug of war between a salt bridge/hydrogen bond (to K68) and halogen bonding (to E114 and V116 backbone oxygens). Here, we test the idea that binding poses might be controllable by pH for ligands with near-neutral pKa, using the conditionally anionic 5,6-DBBt and constitutively anionic TBBt as our models. We characterize the binding by low-volume Differential Scanning Fluorimetry (nanoDSF), Isothermal Calorimetry (ITC), Hydrogen/Deuterium eXchange (HDX), and X-ray crystallography (MX). The data indicate that the ligand pose away from the hinge dominates for the entire tested pH range (5.5-8.5). The insensitivity of the binding mode to pH is attributed to the perturbation of ligand pKa upon binding that keeps it anionic in the ligand binding pocket at all tested pH values. However, a minor population of the ligand, detectable only by HDX, shifts towards the hinge in acidic conditions. Our findings demonstrate that electrostatic (ionic) interactions predominate over halogen bonding.

The Nϵ-Rule for Serine, but Not Cysteine Catalytic Triads.

Catalytic triads, composed of a serine or cysteine nucleophile, a histidine, and a third triad residue (typically Asp/Glu/Asn), are common in enzyme active sites and catalyze a wide variety of chemical reactions. Two types of triads can be distinguished: We refer to them as Nδ- or Nϵ-configured, depending on whether the histidine imidazole Nδ or Nϵ atom is close to the nucleophile Oγ/Sγ. In this study, we have analyzed triad configuration. In structural triads, the more stable Nδ-configuration predominates. For catalytic triads, the configuration depends on the nucleophile. When it is a cysteine residue, both configuration types occur, depending on the family. However, when the nucleophile is a serine residue, the less stable Nϵ-configuration is almost exclusively found. We posit that the energetically less favored conformation is selected for in serine triads to facilitate the otherwise difficult proton transfer from the nucleophile to the histidine residue.

Halogen Atoms in the Protein-Ligand System. Structural and Thermodynamic Studies of the Binding of Bromobenzotriazoles by the Catalytic Subunit of Human Protein Kinase CK2.

Binding of a family of brominated benzotriazoles to the catalytic subunit of human protein kinase CK2 (hCK2α) was used as a model system to assess the contribution of halogen bonding to protein-ligand interaction. CK2 is a constitutively active pleiotropic serine/threonine protein kinase that belongs to the CMGC group of eukaryotic protein kinases (EPKs). Due to the addiction of some cancer cells, CK2 is an attractive and well-characterized drug target. Halogenated benzotriazoles act as ATP-competitive inhibitors with unexpectedly good selectivity for CK2 over other EPKs. We have characterized the interaction of bromobenzotriazoles with hCK2α by X-ray crystallography, low-volume differential scanning fluorimetry, and isothermal titration calorimetry. Properties of free ligands in solution were additionally characterized by volumetric and RT-HPLC measurements. Thermodynamic data indicate that the affinity increases with bromo substitution, with greater contributions from 5- and 6-substituents than 4- and 7-substituents. Except for 4,7-disubstituted compounds, the bromobenzotriazoles adopt a canonical pose with the triazole close to lysine 68, which precludes halogen bonding. More highly substituted benzotriazoles adopt many additional noncanonical poses, presumably driven by a large hydrophobic contribution to binding. Some noncanonical ligand orientations allow the formation of halogen bonds with the hinge region. Consistent with a predominantly hydrophobic interaction, the isobaric heat capacity decreases upon ligand binding, the more so the higher the substitution.

Crystal structures of the EVE-HNH endonuclease VcaM4I in the presence and absence of DNA.

Many modification-dependent restriction endonucleases (MDREs) are fusions of a PUA superfamily modification sensor domain and a nuclease catalytic domain. EVE domains belong to the PUA superfamily, and are present in MDREs in combination with HNH nuclease domains. Here, we present a biochemical characterization of the EVE-HNH endonuclease VcaM4I and crystal structures of the protein alone, with EVE domain bound to either 5mC modified dsDNA or to 5mC/5hmC containing ssDNA. The EVE domain is moderately specific for 5mC/5hmC containing DNA according to EMSA experiments. It flips the modified nucleotide, to accommodate it in a hydrophobic pocket of the enzyme, primarily formed by P24, W82 and Y130 residues. In the crystallized conformation, the EVE domain and linker helix between the two domains block DNA binding to the catalytic domain. Removal of the EVE domain and inter-domain linker, but not of the EVE domain alone converts VcaM4I into a non-specific toxic nuclease. The role of the key residues in the EVE and HNH domains of VcaM4I is confirmed by digestion and restriction assays with the enzyme variants that differ from the wild-type by changes to the base binding pocket or to the catalytic residues.

Protein Domain Guided Screen for Sequence Specific and Phosphorothioate-Dependent Restriction Endonucleases.

Modification dependent restriction endonucleases (MDREs) restrict modified DNA, typically with limited sequence specificity (∼2-4 bp). Here, we focus on MDREs that have an SRA and/or SBD (sulfur binding domain) fused to an HNH endonuclease domain, cleaving cytosine modified or phosphorothioated (PT) DNA. We independently characterized the SBD-SRA-HNH endonuclease ScoMcrA, which preferentially cleaves 5hmC modified DNA. We report five SBD-HNH endonucleases, all recognizing GpsAAC/GpsTTC sequence and cleaving outside with a single nucleotide 3' stagger: EcoWI (N7/N6), Ksp11411I (N5/N4), Bsp305I (N6/N4-5), Mae9806I [N(8-10)/N(8-9)], and Sau43800I [N(8-9)/N(7-8)]. EcoWI and Bsp305I are more specific for PT modified DNA in Mg2+ buffer, and promiscuous with Mn2+. Ksp11411I is more PT specific with Ni2+. EcoWI and Ksp11411I cleave fully- and hemi-PT modified oligos, while Bsp305I cleaves only fully modified ones. EcoWI forms a dimer in solution and cleaves more efficiently in the presence of two modified sites. In addition, we demonstrate that EcoWI PT-dependent activity has biological function: EcoWI expressing cells restrict dnd+ GpsAAC modified plasmid strongly, and GpsGCC DNA weakly. This work establishes a framework for biotechnology applications of PT-dependent restriction endonucleases (PTDRs).

Restriction endonucleases that cleave RNA/DNA heteroduplexes bind dsDNA in A-like conformation.

Restriction endonucleases naturally target DNA duplexes. Systematic screening has identified a small minority of these enzymes that can also cleave RNA/DNA heteroduplexes and that may therefore be useful as tools for RNA biochemistry. We have chosen AvaII (G↓GWCC, where W stands for A or T) as a representative of this group of restriction endonucleases for detailed characterization. Here, we report crystal structures of AvaII alone, in specific complex with partially cleaved dsDNA, and in scanning complex with an RNA/DNA hybrid. The specific complex reveals a novel form of semi-specific dsDNA readout by a hexa-coordinated metal cation, most likely Ca2+ or Mg2+. Substitutions of residues anchoring this non-catalytic metal ion severely impair DNA binding and cleavage. The dsDNA in the AvaII complex is in the A-like form. This creates space for 2'-OH groups to be accommodated without intra-nucleic acid steric conflicts. PD-(D/E)XK restriction endonucleases of known structure that bind their dsDNA targets in the A-like form cluster into structurally similar groups. Most such enzymes, including some not previously studied in this respect, cleave RNA/DNA heteroduplexes. We conclude that A-form dsDNA binding is a good predictor for RNA/DNA cleavage activity.

Crystal structure of the EcoKMcrA N-terminal domain (NEco): recognition of modified cytosine bases without flipping.

EcoKMcrA from Escherichia coli restricts CpG methylated or hydroxymethylated DNA, and may act as a barrier against host DNA. The enzyme consists of a novel N-terminal specificity domain that we term NEco, and a C-terminal catalytic HNH domain. Here, we report that NEco and full-length EcoKMcrA specificities are consistent. NEco affinity to DNA increases more from hemi- to full-methylation than from non- to hemi-methylation, indicating cooperative binding of the methyl groups. We determined the crystal structures of NEco in complex with fully modified DNA containing three variants of the Y5mCGR EcoKMcrA target sequence: C5mCGG, T5mCGA and T5hmCGA. The structures explain the specificity for the two central base pairs and one of the flanking pairs. As predicted based on earlier biochemical experiments, NEco does not flip any DNA bases. The proximal and distal methyl groups are accommodated in separate pockets. Changes to either pocket reduce DNA binding by NEco and restriction by EcoKMcrA, confirming the relevance of the crystallographically observed binding mode in solution.

A protein architecture guided screen for modification dependent restriction endonucleases.

Modification dependent restriction endonucleases (MDREs) often have separate catalytic and modification dependent domains. We systematically looked for previously uncharacterized fusion proteins featuring a PUA or DUF3427 domain and HNH or PD-(D/E)XK catalytic domain. The enzymes were clustered by similarity of their putative modification sensing domains into several groups. The TspA15I (VcaM4I, CmeDI), ScoA3IV (MsiJI, VcaCI) and YenY4I groups, all featuring a PUA superfamily domain, preferentially cleaved DNA containing 5-methylcytosine or 5-hydroxymethylcytosine. ScoA3V, also featuring a PUA superfamily domain, but of a different clade, exhibited 6-methyladenine stimulated nicking activity. With few exceptions, ORFs for PUA-superfamily domain containing endonucleases were not close to DNA methyltransferase ORFs, strongly supporting modification dependent activity of the endonucleases. DUF3427 domain containing fusion proteins had very little or no endonuclease activity, despite the presence of a putative PD-(D/E)XK catalytic domain. However, their expression potently restricted phage T4gt in Escherichia coli cells. In contrast to the ORFs for PUA domain containing endonucleases, the ORFs for DUF3427 fusion proteins were frequently found in defense islands, often also featuring DNA methyltransferases.

Crystal Structure and Directed Evolution of Specificity of NlaIV Restriction Endonuclease.

Specificity engineering is challenging and particularly difficult for enzymes that have the catalytic machinery and specificity determinants in close proximity. Restriction endonucleases have been used as a paradigm for protein engineering, but successful cases are rare. Here, we present the results of a directed evolution approach to the engineering of a dimeric, blunt end cutting restriction enzyme NlaIV (GGN/NCC). Based on the remote similarity to EcoRV endonuclease, regions for random mutagenesis and in vitro evolution were chosen. The obtained variants cleaved target sites with an up to 100-fold kcat/KM preference for AT or TA (GGW/WCC) over GC or CG (GGS/SCC) in the central dinucleotide step, compared to the only ~17-fold preference of the wild-type enzyme. To understand the basis of the increased specificity, we determined the crystal structure of NlaIV. Despite the presence of DNA in the crystallization mix, the enzyme crystallized in the free form. We therefore constructed a computational model of the NlaIV-DNA complex. According to the model, the mutagenesis of the regions that were in the proximity of DNA did not lead to the desired specificity change, which was instead conveyed in an indirect manner by substitutions in the more distant regions.

Crystal structure of the modification-dependent SRA-HNH endonuclease TagI.

TagI belongs to the recently characterized SRA-HNH family of modification-dependent restriction endonucleases (REases) that also includes ScoA3IV (Sco5333) and TbiR51I (Tbis1). Here, we present a crystal structure of dimeric TagI, which exhibits a DNA binding site formed jointly by the nuclease domains, and separate binding sites for modified DNA bases in the two protomers. The nuclease domains have characteristic features of HNH/ßßα-Me REases, and catalyze nicks or double strand breaks, with preference for /RY and RYN/RY sites, respectively. The SRA domains have the canonical fold. Their pockets for the flipped bases are spacious enough to accommodate 5-methylcytosine (5mC) or 5-hydroxymethylcytosine (5hmC), but not glucosyl-5-hydroxymethylcytosine (g5hmC). Such preference is in agreement with the biochemical determination of the TagI modification dependence and the results of phage restriction assays. The ability of TagI to digest plasmids methylated by Dcm (C5mCWGG), M.Fnu4HI (G5mCNGC) or M.HpyCH4IV (A5mCGT) suggests that the SRA domains of the enzyme are tolerant to different sequence contexts of the modified base.

Activity and structure of EcoKMcrA.

Escherichia coli McrA (EcoKMcrA) acts as a methylcytosine and hydroxymethylcytosine dependent restriction endonuclease. We present a biochemical characterization of EcoKMcrA that includes the first demonstration of its endonuclease activity, small angle X-ray scattering (SAXS) data, and a crystal structure of the enzyme in the absence of DNA. Our data indicate that EcoKMcrA dimerizes via the anticipated C-terminal HNH domains, which together form a single DNA binding site. The N-terminal domains are not homologous to SRA domains, do not interact with each other, and have separate DNA binding sites. Electrophoretic mobility shift assay (EMSA) and footprinting experiments suggest that the N-terminal domains can sense the presence and sequence context of modified cytosines. Pyrrolocytosine fluorescence data indicate no base flipping. In vitro, EcoKMcrA DNA endonuclease activity requires Mn2+ ions, is not strictly methyl dependent, and is not observed when active site variants of the enzyme are used. In cells, EcoKMcrA specifically restricts DNA that is modified in the correct sequence context. This activity is impaired by mutations of the nuclease active site, unless the enzyme is highly overexpressed.

Crystal structure of human Acinus RNA recognition motif domain.

Acinus is an abundant nuclear protein involved in apoptosis and splicing. It has been implicated in inducing apoptotic chromatin condensation and DNA fragmentation during programmed cell death. Acinus undergoes activation by proteolytic cleavage that produces a truncated p17 form that comprises only the RNA recognition motif (RRM) domain. We have determined the crystal structure of the human Acinus RRM domain (AcRRM) at 1.65 Å resolution. It shows a classical four-stranded antiparallel ß-sheet fold with two flanking α-helices and an additional, non-classical α-helix at the C-terminus, which harbors the caspase-3 target sequence that is cleaved during Acinus activation. In the structure, the C-terminal α-helix partially occludes the potential ligand binding surface of the ß-sheet and hypothetically shields it from non-sequence specific interactions with RNA. Based on the comparison with other RRM-RNA complex structures, it is likely that the C-terminal α-helix changes its conformation with respect to the RRM core in order to enable RNA binding by Acinus.

The Y. bercovieri Anbu crystal structure sheds light on the evolution of highly (pseudo)symmetric multimers.

Ancestral ß-subunit (Anbu) is homologous to HslV and 20S proteasomes. Based on its phylogenetic distribution and sequence clustering, Anbu has been proposed as the "ancestral" form of proteasomes. Here, we report biochemical data, small-angle X-ray scattering results, negative-stain electron microscopy micrographs and a crystal structure of the Anbu particle from Yersinia bercovieri (YbAnbu). All data are consistent with YbAnbu forming defined 12-14 subunit multimers that differ in shape from both HslV and 20S proteasomes. The crystal structure reveals that YbAnbu subunits form tight dimers, held together in part by the Anbu specific C-terminal helices. These dimers ("protomers") further assemble into a low-rise left-handed staircase. The lock-washer shape of YbAnbu is consistent with the presence of defined multimers, X-ray diffraction data in solution and negative-stain electron microscopy images. The presented structure suggests a possible evolutionary pathway from helical filaments to highly symmetric or pseudosymmetric multimer structures. YbAnbu subunits have the Ntn-hydrolase fold, a putative S1 pocket and conserved candidate catalytic residues Thr1, Asp17 and Lys32(33). Nevertheless, we did not detect any YbAnbu peptidase or amidase activity. However, we could document orthophosphate production from ATP catalyzed by the ATP-grasp protein encoded in the Y. bercovieri Anbu operon.

Halogen bonds involved in binding of halogenated ligands by protein kinases.

Analysis of 664 known structures of protein kinase complexes with halogenated ligands revealed 424 short contacts between a halogen atom and a potential protein X-bond acceptor, the topology and geometry of which were analyzed according to the type of a halogen atom (X = Cl, Br, I) and a putative protein X-bond acceptor. Among 236 identified halogen bonds, the most represented ones are directed to backbone carbonyls of the hinge region and may replace the pattern of ATP-like hydrogen bonds. Some halogen-π interactions with either aromatic residues or peptide bonds, that accompany the interaction with the hinge region, may possibly enhance ligand selectivity. Interestingly, many of these halogen-π interactions are bifurcated. Geometrical preferences identify iodine as the strongest X-bond donor, less so bromine, while virtually no such preferences were observed for chlorine; and a backbone carbonyl as the strongest X-bond acceptor. The presence of a halogen atom in a ligand additionally affects the properties of proximal hydrogen bonds, which according to geometrical parameters get strengthened, when a nitrogen of a halogenated ligand acts as the hydrogen bond donor.

On the role of steric clashes in methylation control of restriction endonuclease activity.

Restriction-modification systems digest non-methylated invading DNA, while protecting host DNA against the endonuclease activity by methylation. It is widely believed that the methylated DNA would not 'fit' into the binding site of the endonuclease in the productive orientation, and thus steric clashes should account for most of the protection. We test this concept statistically by grafting methyl groups in silico onto non-methylated DNA in co-crystal structures with restriction endonucleases. Clash scores are significantly higher for protective than non-protective methylation (P < 0.05% according to the Wilcoxon rank sum test). Structural data alone are sufficient to distinguish between protective and non-protective DNA methylation with 90% confidence and decision thresholds of 1.1 Å and 48 Å(3) for the most severe distance-based and cumulative volume-based clash with the protein, respectively (0.1 Å was deducted from each interatomic distance to allow for coordinate errors). The most severe clashes are more pronounced for protective methyl groups attached to the nitrogen atoms (N6-methyladenines and N4-methylcytosines) than for C5-methyl groups on cytosines. Cumulative clashes are comparable for all three types of protective methylation.

High resolution structure of an M23 peptidase with a substrate analogue.

LytM is a Staphylococcus aureus autolysin and a homologue of the S. simulans lysostaphin. Both enzymes are members of M23 metallopeptidase family (MEROPS) comprising primarily bacterial peptidoglycan hydrolases. LytM occurs naturally in a latent form, but can be activated by cleavage of an inhibitory N-terminal proregion. Here, we present a 1.45 Å crystal structure of LytM catalytic domain with a transition state analogue, tetraglycine phosphinate, bound in the active site. In the electron density, the active site of the peptidase, the phosphinate and the "diglycine" fragment on the P1' side of the transition state analogue are very well defined. The density is much poorer or even absent for the P1 side of the ligand. The structure is consistent with the involvement of His260 and/or His291 in the activation of the water nucleophile and suggests a possible catalytic role for Tyr204, which we confirmed by mutagenesis. Possible mechanisms of catalysis and the structural basis of substrate specificity are discussed based on the structure analysis.

Crystal structure of the antimicrobial peptidase lysostaphin from Staphylococcus simulans.

Staphylococcus simulans biovar staphylolyticus lysostaphin efficiently cleaves Staphylococcus aureus cell walls. The protein is in late clinical trials as a topical anti-staphylococcal agent, and can be used to prevent staphylococcal growth on artificial surfaces. Moreover, the gene has been both stably engineered into and virally delivered to mice or livestock to obtain resistance against staphylococci. Here, we report the first crystal structure of mature lysostaphin and two structures of its isolated catalytic domain at 3.5, 1.78 and 1.26 Å resolution, respectively. The structure of the mature active enzyme confirms its expected organization into catalytic and cell-wall-targeting domains. It also indicates that the domains are mobile with respect to each other because of the presence of a highly flexible peptide linker. The high-resolution structures of the catalytic domain provide details of Zn(2+) coordination and may serve as a starting point for the engineering of lysostaphin variants with improved biotechnological characteristics. STRUCTURED DIGITAL ABSTRACT: lysostaphin by x-ray crystallography (1, 2).

Structural basis of the methylation specificity of R.DpnI.

R.DpnI consists of N-terminal catalytic and C-terminal winged helix domains that are separately specific for the Gm6ATC sequences in Dam-methylated DNA. Here we present a crystal structure of R.DpnI with oligoduplexes bound to the catalytic and winged helix domains and identify the catalytic domain residues that are involved in interactions with the substrate methyl groups. We show that these methyl groups in the Gm6ATC target sequence are positioned very close to each other. We further show that the presence of the two methyl groups requires a deviation from B-DNA conformation to avoid steric conflict. The methylation compatible DNA conformation is complementary with binding sites of both R.DpnI domains. This indirect readout of methylation adds to the specificity mediated by direct favorable interactions with the methyl groups and solvation/desolvation effects. We also present hydrogen/deuterium exchange data that support 'crosstalk' between the two domains in the identification of methylated DNA, which should further enhance R.DpnI methylation specificity.

Crystal structure of the 5hmC specific endonuclease PvuRts1I.

PvuRts1I is a prototype for a larger family of restriction endonucleases that cleave DNA containing 5-hydroxymethylcytosine (5hmC) or 5-glucosylhydroxymethylcytosine (5ghmC), but not 5-methylcytosine (5mC) or cytosine. Here, we report a crystal structure of the enzyme at 2.35 Å resolution. Although the protein has been crystallized in the absence of DNA, the structure is very informative. It shows that PvuRts1I consists of an N-terminal, atypical PD-(D/E)XK catalytic domain and a C-terminal SRA domain that might accommodate a flipped 5hmC or 5ghmC base. Changes to predicted catalytic residues of the PD-(D/E)XK domain or to the putative pocket for a flipped base abolish catalytic activity. Surprisingly, fluorescence changes indicative of base flipping are not observed when PvuRts1I is added to DNA substrates containing pyrrolocytosine in place of 5hmC (5ghmC). Despite this caveat, the structure suggests a model for PvuRts1I activity and presents opportunities for protein engineering to alter the enzyme properties for biotechnological applications.